The incommensurate ground state 
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Double-layer quantum Hall systems possess interlayer 
phase coherence at sufficiently small layer separations, even 
without interlayer tunneling. When interlayer tunneling is 
present, application of a sufficiently strong in-plane magnetic 
field B\\ > Be drives a commensurate-incommensurate (CI) 
transition to an incommensurate soliton-lattice (SL) state. 
We calculate the Hartree-Fock ground-state energy of the SL 
state for all values of By within a gradient approximation, and 
use it to obtain the anisotropic SL stiffness, the Kosterlitz- 
Thouless melting temperature for the SL, and the SL mag- 
netization. The in-plane differential magnetic susceptibility 
diverges as l-By — -Bc|~^ when the CI transition is approached 
from the SL state. 

PACS numbers: 73.40.Hm, 64.70. Rh, 71.10.Pm, 71.45. Gm 



I. INTRODUCTION 

At sufficiently low particle densities and small layer 
separations, double-layer quantum Hall (2LQH) systems 
at total filling factor unitw (vt — 1) can be described 
as broken-symmetry stateil that possess interlayer phase 
coherence, even in the absence of interlayer tunneling.EI 
The 2LQH system can be mapped to an equivalent spin- 
1/2 system by equating "up" ("down'U nseudospins with 
electrons in the upper (lower) layer.l3~El (The electrons 
are assumed to have their real spins polarized.) The ex- 
change energy between the electrons produces a pseu- 
dospin stiffness ps that seeks to keep the pseudospins 
aligned locally. At finite layer separation d, the di- 
rect (Hartree) energy produces a local capacitive charg- 
ing energy that is minimized when the two layers have 
equal electron density. Thus the expectation value of the 
z component of the pseudospin vanishes and the pseu- 
dospin system has an "easy-plane" anisotronv that gives 
the itinerant ferromagnet an XY symmetryQu (in the ab- 
sence of interlayer tunneling). The expectation value of 
a pseudospin at location r can therefore be specified by 
its angle Q(v) in the xy plane. 

In the absence of interlayer tunneling, the 2LQH sys- 
tem picks out a particular global value of B in the ground 
state, spontaneously breaking the (7(1) symmetry of the 
XY ferromagnet. This gives rise to a, linearly dispersing 
Goldstone mode at long wavelengths O'EI Recent measure- 
ments of the zero-bias tunneling conductance in 2LQH 



systems show a huge resonajit enhancement when in- 
terlayer coherence is present This enhancement is re- 
lated to the Goldstone mode of the broken J7(l) symme- 
try, and it has been proposed that the dispersion of the 
Goldstone mode can be observed in tunneling conduc- 
tance measurements by applying a weak parallel mag- 
netic field.Q The XY pseudomagnet also possesses vortex 
excitations called "merons" ; unlike those in an ordinary 
ferromagnet, these vortices are electrically charged, and 
the lowest-energy charged excitationSjJif the system con- 
sist of vorticity-neutral meron pairs .u'El There is also a 
finite-temperature Kosterlitz-Thouless (KT) phase tran- 
sition due to the XY symmetry and the finite pseudospin 
stiffness of the ferromagnet. 

In the presence of interlayer tunneling, the C/(l) in- 
variance associated with conservation of the charge dif- 
ference between the two layers is lost. The finite inter- 
layer tunneling t of the electrons acts as an effective Zee- 
man pseudofield of magnitude 2t along the pseudospin 
X axis and aligns the pseudospins, so that = in the 
ground state. The Goldstone mode disappears, and the 
collective mode becomes gapped. In the presence of inter- 
layer tunneling, merons of opposite vorticity are bound 
together by a "string" that produces a linear attraction 
between the raegons, eliminating the finite-temperature 
KT transitionMS 

Murphy et al. have investigated the effect of an in- 
plane magnetic field on 2LQH systems.tl By varying 
-B|| and studying the energy gap obtained from activation 
energy measurements of the longitudinal resistivity, they 
find evidence for a phase transition between two compet- 
ing QH ground states at a critical value J5|| = Be. These 
two ground states are understood in the pseudospin lan- 
guage as being due to a competition between jdhe tun- 
neling energy t and the exchange energy ps-uB Read 
has presented an appealing analysis of charged (meron 
pair) excitations in this system, focusing on the value of 
the charge-gap near the commensurate-incommensurate 
transition.liy 

The presence of periodically shifts the phase of 
the tunneling matrix elements, resulting in an effective 
Zeeman field for the pseudospins that rotates along the 
planar direction perpendicular to By, with a wavelength 
0o/S[|d, which is the distance required to contain one 
flux quantum (pQ — h/e between the layers. The net 
result is that for gradual variations of the pseudospin 
texture (gradual on the scale of the magnetic length I = 



1 



y/hJeBZ, where B±^ is the strength of the magnetic field 
normal to the plane), the energy of the XY ps«udpspin 
system has the Pokrovsky-Talapov (PT) form,y^t3 
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(1.1) 



up to a constant, where Q = (27rc?/0Q)B|| xz defines the 
parallel magnetic-field wave vector. 



t — tne 



(1.2) 



is the tunneliBg energy (where to is the tunneling energy 
when Q = 0),L3 and 



Ps = (1 - ml)pi 



(1.3) 



is the pseudospin stiffness within the Hartree-Fock Ap- 
proximation (HFA). Here uiz = vi — V2 is the layer im- 
balance, and 
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is the interlayer exchange stiffness when the layers are 
balanced: vi — 1^2 — 1/2, or mz — 0. The layer separa- 
tion is d, and Vj is the filling factor of layer j . 

For small Q (small ), the ground-state energy is min- 
imized by having the pseudospins point in the direction 
of the local (rotating) pseudospin Zeeman field, so that 
9{y) — — Q • r. This is the commensurate ground state, 
and it minimizes the pseudospin Zeeman (tunneling) en- 
ergy. However, it does so at the expense of the exchange 
energy, since neighboring pseudospins are no longer par- 
allel. In the limit of large Q, the cost in exchange en- 
ergy for the pseudospins to align with the rapidly ro- 
tating Zeeman field is prohibitive, and the pseudospins 
become (nearly) uniformly polarized (constant 9), just 
as if i — > 0. The state with uniformly polarized pseu- 
dospin is the large-Q limit of the incommensurate state. 
It turns out that for all finite B\\ > the translational 
symmetry of the pseudospin polarization is broken, and 
a soliton-lattice (SL) state is obtained in the incommen- 
surate phase. 

The SL phase of the PT model can also undergo a sep- 
arate finite-temperature dislocation-mediated KX, tran- 
sition that restores the translational symmetry.t^ This 
work focuses on calculating the ground-state properties 
of the SL state, for all > Be- Interestingly, it is not 
necessary to determine QJj:) .in. order to calculate the to- 
tal energy of the system.tjOEj From the total energy, we 
calculate thermodynamic quantities such as the SL stiff- 
nesses, extending the results of Ref. ^ for the stiffnesses 
and the KT temperature to all i^y. We also calculate 
the SL contribution to the magnetization and suscepti- 
bility, and discuss some possibilities of measuring these 
quantities experimentally. 

The plan of this paper is as follows. In Sec. ||, we 
discuss the single-soliton solutions that follow from the 



equation of motion obtained from the PT energy 
the solitons proliferate; in Sec 
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[II 



For 

the interac- 



tion between soliton lines is discussed, and the separa- 
tion between solitons as a function of Q/Qc is derived. 
In Sec. IV, the compressional and shear elastic constants 
are analyzed, and an estimate is made of the Kosterlitz- 
Thouless temperature for melting the soliton lattice, as a 
function of Q. The interlayer phase coherent 2LQH state 
has a dianoagnetic response to an applied in-plane mag- 
netic field;E3 Section ^ gives a calculation of the in-plane 
magnetization due to the 2LQH state, as a function of 
Q. We conclude with a summary of our results. 



II. SINGLE SOLITONS 



When Q > Qc, i t is convenient to define 9{r) 
Q • r, so that Eq. (1.1) becomes 



Ips {^d-Q]' 
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(2.1) 



Minimizing £ with respect to variations in 9 gives the 
sine-Gordon equation: 



9~ 1 

V^e = -Tsin( 

e 



(2.2) 



where \f2n^p^t. We shall give numerical values for 

our results for a hypothetical "typical" GaAs (effective 
mass TO* K, 0.07me and relative dielectric constant k, 
13) 2LQH sample, which for the sake of definiteness we 
assign the following sample parameters: total particle 
areal density ut = 1.0 x 10^"'^ cm^^, layer (midwell to 
midwell) separation d = 20 nm, and tunneling energy 
to = 0.1 meV. Such a sample would have I w 12.6 nm, 
dji K, 1.6, ~ 6.9 meV for i^t = 1, and e^j^-Kti k, 
8.8 meV. In the HFA, p^ k, 0.03 meV and ^ « 17 nm. 

The commensurate state minimizes the tunneling en- 
ergy by having 9(v) — 0, so that the phase angle 0(r) = 
Q • r follows the tumbling Zeeman pseudofield. The en- 
ergy per area of this state is pj^^ jl — tjl-Kl^ . In the limit 
of large Q, the incommensurate state with constant 9 has 
a lower energy per area, equal to zero. These energies are 
plotted as the solid and dashed lines in Fig. |l|. We there- 
fore estimate that there is a phase transition near the 
point where the commensurate-state and the constant-0 
incommensurate-state energies are equal, at ~ \/2- It 
turns out, however, that the incommensurate state lowers 
its energy by breaking translation invariance, so that at 
finite t and Q > Qc, the value of 9 depends on position. 

Equation (2^) possesses soliton solutions. To see this, 
let us seek solutions of the form 



[ei ■ (r - ro)] 



(2.3) 



wher e ei could be any unit vector in the xy plane. Then 
Eq. (2^) becomes 
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dfe = sin 



(2.4) 



Note that this equation is equivalent to the equation of 
motion for a pendulum of length I, df(f) = —{g/l)sm(j), 
if we replace 9 ^ <j) — tt , ei-r— >t (time), ps — > I, and 
(tunneling amplitude) 1/21:^ g. This analogy is very 
useful in finding the soliton solutions for the PT model. 
In particular, the pendulum can make full circles in a 
given direction. This corresponds to the SL state in the 
2LQH system. 

In analogy with the pendulum problem, we may define 
a conserved quantity analogous to the total (kinetic plus 
potential) energy of a pendulum: 



(1 



COS! 



(2.5) 



Diffe rent iating Eq. (2^) with respect to ei • r and using 
Eq. (2.4) shows that dic = 0, so that c is a constant of 
the motion. Defining (3 = 9/2 then leads to the equation 



di(3 = ±- 



- sin^ 13. 



(2.6) 



It is straightforward to solve Eq. ( |2.6| ) when c = 
by writing / = tan(/3/2), so that dif = ±//^, giving 
9 = 9ss{r), where 



9ss{r) = 4arctan \ e^^^ ^'y^ 



(2.7) 



represents a single soliton in the ei-direction, centered 
at ei • Tq, with width This is shown in Fig. ||. Note 
that 0ss(r) changes by 2tt as ei • r goes from — oo to oo. 
This corresponds to the motion of a pendulum that just 
barely completes a full revolution, and whose period goes 
to infinity. 

The energy per length of a si ngle solito n ma y be com- 
puted by substituting Eq. ( |2.7D into Eq. (2T) and sub- 
tracting the commensurate-state {9 = 0) energy, to ob- 
tain 



^ss / 8 

= Ps ( ^ T 27rei • Q 



(2.8) 



where L2 is the sample length in the planar direction per- 
pendicular to §1. The lowest (soliton) and highest (anti- 
soliton) energy solutions occur for ei = ±Q. Since soli- 
tons in the — Q directions are equivalent to anti-solitons 
in the Q, we shall speak only about solitons with orienta- 
tions defined by ei. The lowest-energy soliton state has 
ei = Q, and its energy per length is 



The value of the critical wave vector Qc will depend on 
the layer imbalance niz = i^i — 1^2- Equations (^_j) and 
(O) give 



Qc(m,) = (l-m2)-i/4Q^(0) 



(2.11) 



in the HFA, where Qc{0) is the value of Qc w hen the 
layers are balanced {niz = 0). Equation (2.11) implies 
that the value of Qc where the CI transition occurs could 
be fine-tuned by adjusting the layer imbalance iriz - i.e., 
by adjusting the gate voltages on the 2LQH sample. Such 
a procedure might be very useful in investigating the CI 
transition. 

When Q < Qc and t > 0, the lowest-energy charged 
excitations are finite-length soliton lines with charged 
meron ends - i.e., charged vortices bound by a soliton 
"string" whose tension is given by Eq. (p!s|). As Q 
increases, the soliton-line "string tension" gets weaker, 
so that the Coulomb repulsion of the merons stretches 
out the string and lowers the energy of the charged 
excitation.El At Q = Qc the soliton-line "string tension" 
goes to zero, and it costs zero energy to make infinitely 
long soliton lines. Since the creation energy per length 
for a soliton decreases linearly with Q (with i3||) for 
Q > Qc, it becomes energetically favorable to form soli- 
tons. The number of solitons created is determined by 
the competition between the (negative) creation energy 
per soliton versus the repulsive (positive energy) interac- 
tions between solitons. Note that Qc < v^/C (the value 
of Q at which the commensurate-state and constant-^ 
incommensurate state energies are equal), so that for 
Q > Qc it is energetically favorable to create solitons. 
Because the solitons are weakly repulsive, the result is a 
soliton-lattice state that we describe below and illustrate 
in Fig. ^. An analogous effect occurs in long Josephson 
junctions, where application of a magnetic field parallel 
to two superconducting planes in close proximity pro- 
duces 2Tr, tari sts in the Josephson phase and generates a 
SL state.llZHia 

To summarize, for Q < Qc, we obtain the commensu- 
rate phase in which 9{r) = 0. For Q = Qc, we introduce 
a single soliton, corresponding to the marginal case of a 
pendulum that makes exactly one full revolution and has 
an infinite period of oscillation. For Q > Qc, we obtain 
a soliton lattice, corresponding to a pendulum making 
complete revolutions in one direction. It has been ar- 
gued that the commensurate to incommensurate soliton- 
lattice (CI) transition at Q = Qc can be identified with 
the phase transition between 2LQH ground states seen by 
Murphy and co-workers ;&B we therefore make the identi- 
fication Qc = 2T:Bcd/ (pQ. 



= 2'n-psiQc - Q), 

J^2 

which goes to zero when Q = Qc, where 



t/Ps 



(2.9) 



(2.10) 



III. SOLITON LATTICE 



We shall now use Eq. (2.6) to determine the SL spacing 
Ls- We do this by noting that over one period of the 



3 



soliton lattice, 9 changes by 27r, so that fi changes by tt. 
We therefore express Lg as 



is 



2 



1 d/3 



(3.1) 



/o v^c2 + sin^ P 

where we have used Eq. ( |2.6| ) and have defined rj 
+ 1, and where 



hi(- 



7r/2 



d/3 



^1 — 772 sin^ /3 



(3.2) 



f) + i(l-r/2) 



ln( 



0-1 



f (1 +1^^ + 

is the complete elliptic integral of the first kind.lli We de- 
fine the SL wave vector Qg = (27r/Ls)ei, so that Eq. (3.1) 
may be written in terms of Qg = |Qs| as 



f]K{i]) 



(3.3) 



Note that 77 ^ 1 corresponds to the CI transition, where 
Qs ^ and Q ^ Qc, whereas 77 ^ corresponds to the 
constant-f? i nco mmen sura te state, where Qs Q ~^ 00. 
From Eqs. (3.3) and (3.2), it follows that 



r 1 - 8e-2-/Q.? 

2 Q, 



Qs/Qc ^ 



1 _ 1 h-Q^ 
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(3.4) 



s/Qc 



00 



In order to determine the value Qg of the SL wave 
vector that minimizes the tot al energy, we express the 
ener gy p er area from Eq. (|]^) as an integral over /? [cf. 
Eq. (3.1)] and obtain 



L1L2 e 



- Q • Qs ) e 



(3.5) 



QcQs^ 
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ps [\Q^ - Q • Qs 



< 



QcQs 



1 + 4e ^ 



\ (Q - Qs)' - 1/ (2Qse' 



7/^1 

77^ 



where 





E{^) EE / d/3y 1 - 7;2 sin^ /3 



(3.6) 



1 + 1(1-77^) 



ln(- 



7; ^ 



is the complete elliptic integral of the second kind,0 and 
L\ is the sample length in the planar direction paral- 
lel to §1. Agreement between Eqs. (2.8) and (3.5) in 



the thermodynamic limit is obtained by equating Qg = 
(27r/ii)ei, so that — L\ ■ak, Q ~ Qc- 
All the terms in the 7; — ^ 1 limit (near the CI transition. 



where the solitons are well separated) of Eq. (3.5) have 
simple physical interpretations. The first and last terms, 
Ps((3^/2 — 1/C'), constitute the (commensurate-phase) 
energy per area of the PT model in the absence of soli- 
tons. The second and third terms, PsiQcQs — Q • Qs) are 
just the creation energy per area for the solitons in terms 
of interacting soliton lines [see Eq. (p.9|)]. The fourth 
term, 4ps(3cQs exp(— 27r/Qs'?): is the exponentially weak 
repulsive interaction energy per area between the soli- 
tons. Because Qs/'^tt = 1/Ls is the density of soliton 
lines, the fourth term says that near the CI transiton, 
the interaction energy per length between two parallel, 
straight, and infinitely long solitons lines separated by a 
distance is 



lim — 

Ly 



SnpsQce 



(3.7) 



Hence for Q > Qc soliton lines proliferate rapidly un- 
til the repulsion between the solitons become signifi- 
cant. The notion of an effective repulsive interaction be- 
tween sipe- Gordon solitons was discussed by Perring and 
Skyrme,ta who obtained the exponentially weak repul- 
sion found above. The arguments of Ref. |l^ imply that 
when the solitons are close together (Ls/C ~^ 0) at large 
Q/Qc, the repulsive potential energy per length between 
soliton lines is {tt^ /2)psQc{^/ Lg). This latter repulsion 
is due to boundary condition that 9 must change by 2tt 
over the small distance Lg, which implies a large gradient 
energy. 

Th e va lue of Qs which minimizes the energy per area 
[Eq. ( |3.5[) ] is found by setting to zero 



d 



dQs \L1L2 



Q 



Eiv) 



Qs-Q, 



where we have used the identit; 

dE _ E{ri) - Ki-q) 
dr/ rj 



(3.8) 



(3.9) 



It is not difficult to show that the second derivative of 
the energy per area with respect to Qs is always positive; 
thus the extremum calculated above is a minimum. It 



follows from Eq. (3.8) that the energy is minimized for 
Qs = Q, and for rj — fj such that 



Q _ Ejfj) 

Qc fj 

where Q = |Q|. 

We define the CI closeness parameter 



(3.10) 
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(3.11) 



e=Q/Qe-l, 



which vanishes at the CI transition; from Eqs. (|3.6| ) and 
( |3.10D , it foUows that for 77 -> 1, 



1 — ?7 



In 



1 — 77 



(3.12) 



so that asymptotically, 
' l-2e/ln(l/e), 

\2 Q J ^ i\2 Q J 



s/Qc 



s/Qc — > 00 



(3.13) 



Equations (3^) and ( 3.1C ) together allow us to deter- 
mine the equilibrium SL wave vector Qs(Q) that mini- 
miafis the energy, in terms of the parallel-field wave vector 
QEI We ha ve pl otted this in Fig. |. From Eqs. (Q, 



( pl0| ), and ( pl3| ), it follows that 



Q 



(V2)^ 
K{fl)E{fj) 
r (7rV2)/ln(lA), 



(3.14) 



1 - 



32 Q 



Q/Qc 

Q/Qc 



where the Q Qc limit is true asymptotically.EilLUtill We. 
note however that, as found by Pokrovsky and Talapovtll 
and discussed in Ref. |2^, the meandering of soliton lines 
at finite temperature renormalizes the dependence of the 
soliton-line density on the parallel magnetic field, so that 
Qs oc ^/e sufiiciently close to the CI transition. 

The minimized value £ / L1L2 of t he e nergy per a rea 
at Qs — Qs is found by using Eqs. ( |3.5| ) and ( 3.10 ) in 
Eq. (3.6) to obtain 



1 



t 



L1L2 
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(3.15) 



-ipsQ?(7rV8-l), Q/Qc 

The SL state energy per area is represented by the dots 
in Fig. |l|. 

Although it is not needed for calculating the stiffnesses 
or suscep±|ibility of the SL, the SL solution for ^(r) is 
given byta 

sin[i(^ - T,)] = sn(ei • (r - ro)/???, v), (3.16) 

where sn denotes the sine-pnplitude Jacobian elliptic 
function with parameter ?7.E3 Near the CI transition, 
when Q — > Qc, is very close to being a periodic su- 
perposition of single-soliton solutions, spaced apart by 
Ls, so that ^(r) « ^ss(r — jLs), where 0ss(r) denotes 



the single-soliton solution, Eq. (2.7). In the SL state, 
especially away from the CI transition, it is sometimes 
useful to work with 

0s(r) = ~e{Y) - Qs • r = e{v) + (Q - Qs) • r, (3.17) 

because it is periodic in the SL spacing, so that 6's(r + 
Ls) — 0s(r), where Ls = Qs27r/Qs. In the limit Q/Qc 
oOj Qs Q and ^s ^ — > 0, so that we may regard 
0s (r) as a sm all q uantity. Expressing the sine-Gordon 
equation [Eq. (2^)] in terms of ^s and working to lowest 
order in ^s gives V^^s ~ (1/?^) sin[Qs • (r — Tq)], so that 



lim 0s (r) 

Q— »oo 



sin[Qs-(r-ro)]. (3.18) 



IV. STIFFNESSES OF THE SOLITON LATTICE 



The elastic constants of the soliton lattice are given by 
the stiffness tensor Kij. The change in the energy due 
to varying the spacing between the soliton lines is char- 
acterized by the compressional stiffness Kn. The shear 
stiffness K22 is associated with the change in energy due 
to shearing the upper and lower ends of the soliton lines 
in opposite directions, and is equivalent to a combined ro- 
tation and compression. We use the calculated stiffnesses, 
to describe the _B|| -dependence of the KT temperatureEHl 
for the dislocation-mediated KT melting transitionEJ of 
the soliton-lattice. 

We calculate the stiffness tensor by two methods. 
First, we calculate the stiffness Kij that is obtained by 
differentiating £/LiL2 in Eq. ( ^j| ) with respect to the 
I, j components of Qs for fixed Q at the extremal, where 
Eq. (|3.8|) is zero. Then we calculate the stiffness tensor 



Kij obtained by calculating the effects of fluctuations of 
the angle variable 9{y) away from its ground-state value, 
Eq. ( CT - 

We begin by calculating the stiffness tensor Kij f rom 
the dependence of the ground-state energy [Eq. (3.5)] on 
the soliton-lattice wave vector Qs. The expressions we 
obtain for Kij by this method have been carried out for 
all values of Q > Qc and agree with the results obtained 
in Ref. |l^, in the limit Q ^ Qc- Because the stiffnesses 
involve the second derivative of £/LiL2 with respect to 
the components of Qs at fixed Q, the terms in Eq. (3.5) 
that depend explicitly on Q (including the term — PsQ • 
Qs that gives the orientational dependence of the energy 
per area) do not contribute to Kij. Thus 
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Kij = lim 

Q,^Q, aQsiOQs- 
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= psQc lim 
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Qi 9Qs 



(4.1) 



(4.2) 
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0, 



Q 



where we have used Eqs. (3.8) and (3.10). Since Qs2 = 0, 
it foUows that K12 = K21 



0. 



Using the results of Sec. [II, and the identit 



dK 
drj 



E{v) K{r,) 



ri{l - jf) 
from which it follows that 
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(7r/2)2 (1 - ^2) /Q 
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we find that the compressional elastic constant Kn is 
equal to 



Ell 

Ps 
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Q/Qc 
Q/Qc 



00 



In the limit Q/Qc 1 when the soliton lines are far 
apart, Kn is very small [of order e ~ exp(— Lg/Oi see 
Eq. ( ^.14| )]. The energy cost of compression very close to 
the CI transition is due to the exponentially weak inter- 
soliton interaction energy. The energy per area due to 
the string tension of the soliton lines (the term PsQcQs) 
does not contribute to Kn, although it does contribute 
to K22, as we explain below. 

As explained in Ref. 23, soliton lines meander at fi- 
nite temperature and are no longer straight; collisions 
between meandering soliton lines produce an effective 
entropic repulsion between the solitons that dominates 
the exponential repulsion at any nonzero temperature, 
for Ls/S, sufficiently large. This effect renormalizes the 
compressional stiffness Kn upwards so that it becomes 
proportional to T^. In the limit Q/Qc — > 00, the tunnel- 
ing term in the PT energy (|2.l|) fluctuates on very short 
length scale and averages to zero, so that Eq. (2.1) be- 
comes the isotropic XY model; thus one expects Kn to 
approach the p seud ospin stiffn ess ps in that limit. 

From Eqs. (3.10) and (4.1), it follows that the shear 
elastic constant K22 is given by 



K22 

Ps 



Q_ 

Qs 



(2/7r2)ln(l/e), Q/Qc 



1 



1 



2(40") ' Q/Qc 



(4.6) 



As expected, the shear stiffness K22 approaches the pseu- 
dospin stiffness ps in the limit Q/Qc — > 00. But in the 
limit Q ^ Qc, K22 diverges as Qc/Qs- The origin of 
this effect is that the shear motion described by K22 in- 
creases the total length of the soliton lines. By definition, 
K22 describes a shear displacement in which Qs2 changes, 
while Qsi remains fixed: i.e., the solitons lines tilt away 
from their equilibrium "vertical" (§2) direction by a small 



angle (/) = Qs2/Qs, while keeping their "horizontal sepa- 
ration" ei • Lg = Lg constant. This shear motion re- 
duces the mean soliton separation (the separation along 
the direction perpendicular to the tilted soliton lines) to 
Ls = Ls cos (p, so that Qs Qs/ cos <j). The if22-shear 
corresponds to a global rotation plus a compression of 
the SL. Packing the solitons lines more closely together 
in a fixed sample area produces more total soliton line 
length, which costs more soliton-line creation energy. 

Because the term — p^Q • Qs that contains the orien- 
tational dependence of the energy per area [Eq. (^)] is 
linearly proportional to Qsi, it cannot contribute to the 
stiffnesses (which are proportional to second derivatives 
of the energy per area with respect to the components 
of Qs) at all. It might be supposed that the rotation of 
the soliton lines that occurs upon shearing should cost 
energy, but this is not so, because a shear is a combi- 
nation of rotation and compression, rather than a pure 
rotation. The creation energy per area of the SL near 
the CI transition {Q « Qc) is PsiQcQs — Q • Qs), and 
consists of two terms. The first term (psQcQs) is just the 
total line length of the solitons times the line tension, 
divided by the total area. The second term (— psQ • Qs) 
depends explictly on the angle 4> between Q and Qs, and 
determines the orientation of the SL because it is min- 
imized by choosing Qs along Q (i.e., = 0); thus, a 
different choice of SL orientation (i.e., a pure rotation of 
the soliton lines) would cost more energy. Interestingly, 
the second term in the creation energy is unchanged by 
a shear, because the energy cost of rotating the soliton 
lines is exactly offset by the increase in total soliton line 
length: —[Q{Qs/ coscp)] coscp — ~QQs, independent of 0. 
The only change in the creation energy comes from the 
first term, which depends only on the density of soliton 
lines: PsQdQs/ cos (j)). Sufficiently close to the CI transi- 
tion (i.e., when Lg/^ 3> 1), the exponentially small inter- 
actions may be neglected in comparison to the creation 
energy. For Q — > Qc and small shear {Qs2/Qsi 1), 



L1L2 



PsQcQs - t/2TT£'^ 
1 Qs2 

" 2 Qs 



(4.7) 



t/2Tif-, 



so that i^22 — > PsQc/Qs ^ 00 as Q ^ Qc, in agreement 
with the results of Ref. |l^. 

The fact that bilayer phase-coherent 2LQH states can 
exhibit a finite-temperature Kosterlitz-Thouless (KT) 
transition in the absence of interlayer riunneling {t — > 
0) has been discussed in earlier worka^El Although fi- 
nite t removes the possibility of a KT transition in the 
commensurate phase of 2LQH systems by altering the 
nature of the long-range, interaction between vortices 
("merons" in this case)j3 the SL phase does support 
a finite-temperature KT tra|Bsition due to dislocation- 
mediated melting of the SL.E3 As discussed in Ref. |l^, 
the KT temperature may be estimated as 
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(7r/2)ps Ps ^ Ps 



(4.8) 



4^1^ / (2/7r2)y'eln3(l/e), Q/Qe ^ 1 



(Qse)^'?^ 



1 + i ft) - Q/Qc^oo 



where we have used the zero-temperature values for Kij 
that we calculated previously to make a rough esti- 
mate the KT temperature. As mentioned earlier, finite- 
temperature fluctuation effects can strongly renormalize 
Kiif3 and may also effect K22- Our results agree with 
those of Ref. |lO| in the limit Q ^ Qc- We have plotted 
the compressionai (i^ii) and shear {K22) stiffnesses in 
Fig. H, together with the KT temperature. 

The KT transition would be most easily measured from 
the temperature dependence of the linear response to 
oppositely directed currents in each layer. This would 
require double-layer electron devices with layers that 
could be contacted separately. Unfortunately, the leak- 
age currents produced when the interlayer tunneling is 
not vanishingly small would make it difficult, perhaps 
impossible, to set up oppositely directed currents in each 
layer. However, because the SL dislocations are electri- 
cally charged, it might be possible that the KT transition 
could be signalled by an increase in the usual longitudinal 
resistivity pxx{T), measured in devices with the current 
flowing in the same direction in both layers. The in- 
crease in (T) would originate from the proliferation of 
unbound charged dislocations above the transition tem- 
perature. 

We now calculate an alternate stiffness tensor Kij 
by examining the effect of deviations of the angle vari- 
able 9{r) from its ground-state value. We write 9{r) = 
6*0 (r) -I- 69 (r), where 9o{r) is th e gr ound-state solution 
that minimizes the PT energy (2T) and therefore sat- 
isfies Eq. (2.2), and 66{r) is the deviation of 9 from its 
ground-state value. We do not include dynamics here, 
because our focus is on ground-state, rather than excited- 
state, properties. The PT energy for 9 is i?pT[^o + S9] = 
Ept[9o] + SH, where Ept, given by Eq. (O), is the PT 



energy from which the ground-state 9o{r) is determined 
via Eq. dU), and 
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~ 2 



27r^2 



27r£2 



t cos 9r) 



tcosl 



f 27r^Vs (yS9] 



(4.9) 



where we have kept terms up to quadratic order in 59. 
There are no terms linear in 69 because 9o{r) is deter- 
mined by minimizing EpT with respect to variations in 
9. The total energy is minimized by choosing 69 from 
among the eigen valu es of the bracketed Schrodinger-like 
operator in Eq. (4.9), so that 



t cos 6*0 - 27rrpsV' 



(4.10) 



If we take B|| = _B[|y so that Q = QSc, then 9o{r) de- 
pends only on x, 6H is translationally invariant in the y 
direction, and we may write 69q{r) = e-xjp{iqyy)69q^{x). 
The term icos6'o in Eq. ( 4.10| ) is periodic in the x direc- 
tion, with a period of Z^ . As shown in Ref. |l^, when 
^o(r) has the form ( ^.16 ), Eq. (4.10) becomes Lame's 
equation, after a simple rescaling of x. Lame's equation 
has three simple solutions, two of which have low-energy 
limits. The first type of solution has zero energy and 
corresponds to a uniform translation of the vortex lines, 
69 cx d9Q/dxQ, where Xq is the x component of Tq in 
Eq. ( 3.161 ). Of greatest interest to us are the type of so- 
lutions which have been called "vortex oscillations" in 
the context|-e|f long Josephson junctions in parallel mag- 
netic fields.lla From Ref. ^ it follows that in the long- 
wavelength limit. 



6£ 
L1L2 



- {Kiiql+psql) 



(4.11) 



wher e K n is equal to the compressionai stiffness in 
Eq. (4.5), and q^ is the crystal momentum along the x 
direction. 

Although Kii — Kii, we find that K22 — Ps ^ K22- 
The reason for the discrepancy between K22 and K22 is 
not obvious. It may be that using K22 — Q/Qs for the 
transverse stiffness is valid only at very long wavelengths 
qy < 1/YsLf^ where IsL = Ls\/K^2/K^ is Read's es- 
timate of the minimum u^distance for SL dislocations to 
interact logarithmically.t3 It is also possible that K22, 
which is based on a calculation that assumes 69 is ev- 
erywhere small, is not able to describe uniform shear, 
which can move solitons lines far from their equilbrium 
positions on the scale of the soliton line thickness ^. The 
relationship between K22 and K22 requires further clari- 
fication, including a stronger argument for preferring K22 
over K22 in estimating the KT temperature. 



V. IN-PLANE MAGNETIZATION 

The interlayer phase coherent 2LQII state exhibits an 
in-plane magnetization M|| in the presence of an in-plane 
magnetic field By. The in-plane magnetization can be 
calculating by differentiating the minimized ground-state 
energy per volume with respect to the parallel magnetic 
field: 



Mil 



1 



d£ 



LiL2ddB 



£ 



27r. d 

2 X 

(t>o dQ \L1L2 



(5.1) 



where £ / L1L2 is given by Eq. (3T5). In order to carry 
out t he d ifferentiation in Eq. (5.1), we first differen tiate 
Eq. ( |3.1C| ) with respect to Q and make use of Eq. (3.3) 
to obtain 



dri 
dQ 



-3 Qs 



(5.2) 
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We note here that the Q^£^/4 term in Eq. ( |5.2D arises 
from differentiating Qc in Eq. ( |3.10 ) with respect to Q\ 
the Q-dependence of Qc is due to the dependence of the 
tunnehng matrix element on the parallel magnetic ficLd, 
t = to exp(— Q^£^/4), which is a single-particle effcct.tll 
The tunneling part of the equilibrium energy per area 
will also give a contribution to the in-plane magnetiza- 
tion proportional to dt/dQ, again due to dependence of 
t on Q. It is convenient to separate My into two parts: 
Mj| = MsL + Mt, where Msl is calculated at fixed t (t 
independent of Q), and Mt arises from the Q-dependence 
of i. 



Using Eqs. (5.1) and (5.2) it is straightforward to show 
that the SL contribution to the parallel magnetization is 



where 



MsL = -Mozx(Q-Q,)/Qc, 



Mo = 27rpsQc/0o ^ 0.5A/m 



(5.3) 



(5.4) 



sets the scale of the SL magnetization, and the numerical 
estimate of Mo isgiven for the "typical" GaAs sample 
described in Sec. O. For such a sample, 



Mod = 27rpsQcd/0o ^ 10^^ A 



(5.5) 



sets the scale of the SL magnetic moment per unit area. 
Thus magnetometers with sensitivities in the range of 
10~^^ Am^ to 10~^^ Am^ would require sample areas in 
the range of 1 mm^ and 1 cm^ to measure MgL. The 
magnitude of the SL magnetization behaves like 





-1 


Mo 





l-(7rV2)/ln(l/e), Q/Qc 

64 (■2^ ) ' Q/Qc 



1 



(5.6) 



The SL magnetization may also be calculated directly 
from the pseusdospin supercurrent density,u J^^ = Ji — 
J2 = {2ps/h)V6, and the definition of the magnetic mo- 
ment. The electrical current I in layers 1 and 2 is 



I = Ii = -I2 = Lyieps/h)iV9 - Q) 



(5.7) 



The maspetization produced by the above current is 
thereforeO 



Msl = 



1 



-zx 



LxLyd 
—zx — 



-dxdz 
c 



(5.8) 



(ve - ci)dx 

zx(Q - Qs), 



eft 



in agreement with Eq. (|5^ 

The 2LQH interlayer phase coherent state mayJpe re- 
garded as a pseudospin-channel superconductor^ The 
magnetization M|| is due to pseudospin supercurrents, 
corresponding to electrical currents traveling in oppo- 
site directions in each layer, which partially screen B||. 



For Q < Qc, B|| is maximally excluded from the region 
between the planes; but when Q > Qc, additional 
penetrates the region between the plates in the form of 
solitons that each carry a single flux quantum (f>Q — h/e 
(corresponding to a phase change A9 — 2tt), leading to 
a precipitous decline in the magnetization. The direc- 
tion of M|[ is opposite to B||, in accord with Lenz's law. 
An exactly analogous effect occurs for magnetic fields 
applied parallel to narrow insulating regions (Josephson 
junctions) between superconductors. 

The contribution to the magnetization due to the Q- 
dependence of the tunneling matrix element t is 



(5.9) 



4_Q 

TT Qc 



Qc 



In the commensurate phase, 

Mt = -(A/o/2)(£/^)2zxQ/(5c, and the magnitude of the 
tunneling contribution to the parallel magnetization be- 
haves like 



Mt 



Mn 



l-4/ln(l/e), Q/Qc-^l 
i(7r/4)2(g,/Q), Q/Q, ^00 



(5.10) 

The SL magnetization is plotted in Fig. ^. It is useful 
to compare the total SL magnetization Mq to the scale 
of the Landau diamagnetism in a, vt — ^ 2LQH system: 



Mo 



MQ(t)od 



16^^ ^0.1, 



(5.11) 



where /Xg = eh/2m* is the effective Bohr magneton, and 
we have made use of the parameters for the "typical" 
GaAs sample described in Sec. ||. This shows that the 
SL magnetization is expected to be roughly an order of 
magnitude smaller than the Landau diamagnetism. It 
is interesting to note that the weak signals associated 
with orbital diamagnetism in two-dimensional electron 
systems have been measured at high magnetic fields by 
torsional magnetometry,Ej SQUID magnetpmetryjH3 and 
micromechanical cantilever magnetometry.E3 The torque 
on the 2LQH sample in the presence of both a perpen- 
dicular and parallel magnetic field has both a Landau- 
diamagnetic component t±_ — {LxLyd)M±B\\ and a SL 
component, r|| = {LxLyd)M\\B±. The smallness of the 
parallel moment My is offset by the large perpendicu- 
lar magnetic field B^ in the expression for the torque t\\ . 
The SL magnetization MgL might also be measured using 
high-field magnetometry techniques mentioned above. 

It follows from Eq. (5.1) that when M|| is plotted 
against B\\ , the area under under the resulting curve from 
B\\ = to By = 00 is 



^^^"'^^i' = TTd 

^x^yU' 



[£{B^\=oo)~£{B\\=Q)\ (5.12) 
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\2TTpd 



B„ =0 



where {nix) is the ground-state expectation value of the 
X component of the pseudospin order parameter, which 
has a value of ^1 — in the HFA, so that {nix) = 1 
(in the HFA) when the layers are balanced (m^ = 0). 
Thus the area under the My versus i3|| curve may be re- 
garded either as a measurement of t (if quantum fluctua- 
tions in the ground state are neglected) , or as a measure 
of order-paijameter suppression (of nix) due to quantum 
fluctuationsEa (if t can be measu red separately). 

Equations (5.3) and (3.14) show that Msl = 
— A/qzxQ/Qc — — (xo/mo)B|| in the commensurate 
phase (Q < Qc), where xo = Mo(27r/0o)^Ps'i sets the 
scale of the SL contribution to the magnetic susceptibil- 
ity. The SL contribution to the parallel-field magnetic 
susceptibility is defined for fixed t (independent of Q): 



XSL = Mo 



fdMsL\ ^ 2TTd fdMsL 
[ )^ <Po [ dQ 

1, Q/Qc < 1 

^(7rV2)/[eln(l/e)], Q/Qc - 1+ 



(5.13) 



where we have used Eqs. (3.1C) and (|4 
Xo = Mo(27r/(/)o)^Psd 



(5.14) 



4:TTa tr 



3 X 10"^ 



and the numerical estimate of xo is given for the "typi- 
cal" GaAs sample described in Sec. Here a « 1/137 is 
the fine-structure constant. We have plotted the suscep- 
tibility in Fig. 1^. Note that near the CI transition, the 
susceptibility diverges like 1/(Q — with logarithmic 
corrections. 

It might be possible to measure the SL magnetization 
or even the SL magnetic susceptibility by varying the 
gate voltages of the sample, in order to adjust Qc and 
tune close to the CI transition. Measuring xsl might 
be possible using AC modulation of the gate voltages 
in order to AC modulate the layer imbalance and 
therefore the critical wave vector Q^- By such a method, 
the ratio Q/Qc could be AC modulated just above and 
below unity, allowing xsl to be determined at or near 
the CI transition. As an example, we compute here that 
part of dMsh/diTT-z which is proportional to xsl: 



dnriz 



diUz 
1 



dQc dniz 
4>o dQc Q 



dQc Q OMsL 
Qc dQ 



(5.15) 



fio STrd dmz Qc ''^^^ 
ni^ Mo Q f Ps 



I -ml 2 Qc 



where we have used Eq (2.11) in the last line. Equation 
( |3.15 ) shows that dMsh/dniz has a contribution propor- 
tional to XSL, which diver ges li ke ps/Kn near the CI 
transition. However, Eq. ( 5.15 ) vanishes for balanced 
layers (m^ — 0); thus, layer imbalance is required to 
measure dMsh/dniz- It turns out that when the layers 
are not balanced (niz ^ 0), the solitons have a nonzero 
electric dipole moment per length;Lj this changes the in- 
tersoliton repulsion from being exponentially weak to a 
power law. The net result is that the stiffness K\\ near 
the CI transition is strengthened from being essentially 
linear in {Q — Qc) to being proportional to ^/Q — Qc, so 
that the divergence in xsl becomes an inverse square- 
root singularity. In real samples, we expect this singu- 
larity to be smoothed out by finite temperature and dis- 
order. 

In practice, is varied by adjusting gate voltages. 
For a small change SVg in a gate voltage away from bal- 
ance {mz = 0), the change Sniz in the layer imbalance 
is linearly proportional to SVg- For the "typical" GaAs 
sample described in Sec. |l[ a rough estimate indicates 
that 



SrUz 



2SVa 



cDgut I [treo) 



3 ^Vc/volt, 



(5.16) 



if we take Dg = 10^^ m to be the distance between 
the gate and the double layers. There is an additional 
complication in changing the layer imbalance by adjust- 
ing one of the gate voltages. Unless both the back and 
front gates are adjusted together in a coordinated way, 
the change in gate voltage 5Vg will also change the total 
filling factor ut by an amount 5vt, where 



5vt 



5Vg 



eDGnr /{ereo) 



1.5 6Vg /yolt, (5.17) 



if Dg = 10-^ m. 



VI. SUMMARY 



In the presence of a sufficiently strong parallel mag- 
netic field, a i^T = 1 2LQH device undergoes a transition 
to a soliton-lattice (SL) state. We have investigated the 
ground-state properties of the SL state for all values of 
the parallel magnetic field, with an eye towards possible 
experimentally measurable effects. 

We found that the SL contribution to the orbital mag- 
netization rises in the commensurate phase {Q < Qc) 
with Q, and quickly drops to zero in the incommensu- 
rate phase {Q > Qc). An estimate of the size of the SL 
magnetization shows that it could j-bje detected by sen- 
sitive magnetometry techniques E3c3 The SL magnetic 
susceptibility shows a singularity at Q = Qc, and it was 
proposed that this signature of the CI transition might 
be detected by varying the gate voltages so as tune close 
to the CI transition. 
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The longitudinal and transverse SL stiffnesses were 
computed and used to estimate the temperature of the 
KT transition, which might be indicated experimentally 
by an increase in pxx at the transition. A more sensitive 
signal of the KT transition would be obtained by mea- 
suring the transresistivity (Coulomb drag) as a function 
of temperature, although this would require separately 
contactable layers. However, the leakage currents pro- 
duced by any sizeable interlayer tunneling might make it 
impractical to set up oppositely directed currents in each 
layer. 

In this paper, we have largely neglectejd the existence 
of disorder. As pointed out by FisherEll and discussed 
by Readta, randomness in the tunneling t due, for exam- 
ple, to small variations in the barrier thickness, pins the 
SL domain walls randomly and destroys the long-range 
order, no matter how weak the randomness. This puts 
limits on how closely Q can approach Qc ~ i-C, how small 
e can be. 

We have not discussed here the dynamics-of the soli- 
tons - i.e., the motion ofjindividual solitons£3 or the col- 
lective motion of the SLe3 - which also produce experi- 
mental signatures of the incommensurate SL state. We 
have focused instead on ground-state properties. One of 
us has recently found that when the layer densities are 
made unequal by adjusting the gate voltages of the device 
(keeping the total filling factor equal to one), the layer 
densities in the SL state, become "rippled" , resulting in 
a dipole density wave.EEl Preliminary calculations show 
that the sudden onset of such a rippled SL state may 
give large contributions to the differential capacitance of 
2LQH systems, especisJiy for the interlayer capacitance 
("Eisenstein ratio" ).a'E3 Further work on the properties 
of the rippled state is in progress, but it may be that 
sensitive measurements of differential capacitance in un- 
balanced {vi ^ 1/2) 2LQH systems could detect the CI 
transition. 
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FIG. 2. At Q = Qc, the system admits a single soliton, in 
which ^(r) twists by 2n over a distance 5. 



5.0 



4.0 



3.0 



2.0 



1.0 



0.0 



-100.0 



-50.0 



0.0 



50.0 



100.0 



FIG. 3. Soliton-lattice for Q slightly larger than Qc 



FIG. 1. Energy of the commensurate (C), constant-^ in- 
commensurate(I), and incommensurate SL (S) phases, versus 
the parallel magnetic-field wave vector Q. For Q > Qc, the 
SL phase (dotted) has the lowest energy. 
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FIG. 5. Soliton-lattice stiffnesses px = Ku and py = K22, 
and their geometric mean, which is proportional to the 
Kosterlitz-Thouless temperature, Tkt- 



300.0 f 




FIG. 7. SL contribution to the magnetic susceptibility, 
which diverges as 1/(Q - Qc) ln[Qc/{Q — Qc)] for Q Qc- 



12 



